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Abstract 

We calculate the Compton and annihilation production of a soft 
static lepton pair in a quark-gluon plasma in the two-loop approxima- 
tion. We work in the context of the effective perturbative expansion 
based on the resummation of hard thermal loops. Double counting 
is avoided by subtracting appropriate counterterms. It is found that 
the two-loop diagrams give contributions of the same order as the one- 
loop diagram. Furthermore, these contributions are necessary to ob- 
tain agreement with the naive perturbative expansion in the limit of 
vanishing thermal masses. 
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1 Introduction 



We consider the production of a soft static lepton pair in a quark gluon 
plasma. In our approach we follow strictly the hard thermal loop (HTL) 
scheme of [||, |2|, appropriate for a plasma in equilibrium at high tempera- 
ture with a small coupling constant between the quarks and the gluons in 
the plasma. We thus construct the loop expansion using effective vertices 
and propagators instead of bare ones. The production rate of static virtual 
photons has already been evaluated at the one-loop level in the effective 
theory Q . Previous works [Q, |5| have shown that two-loop diagrams are of 
equal importance as the one-loop diagram. These large contributions are 
associated with new processes (namely bremsstrahlung) which arise only at 
the two-loop level. However, these processes are only a part of the full 
two- loop diagrams. In this paper we extend our previous work || and com- 
plete the calculation of these diagrams. The underlying physical processes 
to be considered are Compton and quark-antiquark annihilation which are 
already present at one loop. The underlying physical processes indicates 
that there is a possible double counting between one and two- loop results. 
More precisely, the two-loop diagrams of Fig. 1 are already contained in the 
one-loop diagram of ||] (with effective propagators and vertices) when the 
gluon is hard and time-like. We take care of this problem by subtracting 
the appropriate counterterms. 

We show that the two-loop contribution corrects in a crucial way the 
calculation of the Compton and annihilation processes based on the one- 
loop approximation || f?j. Adding one and two loop contributions we find 
that in the limit of vanishing thermal masses we recover the result already 
found in ^ where the authors used the bare theory. This solves a long 
standing puzzle related to the fact that the one-loop approximation in the 
effective theory did not reduce to known results in the naive perturbative 
expansion ||. 

It is not possible to derive analytically the full leading expression, in 
powers of g, of the rate of static virtual photon production up to two loops 
in the effective theory. However, it is relatively easy to obtain analytically 
the large logarithmic ln(l/g) terms. After some general considerations con- 
cerning our approach and the approximations useful to extract the leading 
logarithmic behavior we review the one-loop results and then discuss in 
Sec. |4] the two-loop calculation. The next section is devoted to an approxi- 
mate evaluation of the counterterms. Combining everything we then obtain 
in Sec. ^ the rate of virtual photon production. In the Appendix we give the 
exact result for the counterterms, using effective vertices, where we show 
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that the leading logarithmic behavior is the same as that obtained using the 
simplified version. 



2 General considerations 

It is well known that the production rate of a photon of invariant mass 
\/Q 2 decaying into a lepton pair is proportional to the imaginary part of 
the retarded vacuum polarization of the considered photon [QQ, nU: 



dN 
dtd 3 x 



dqod?q a 



n s (<2o)ImII /(<?o,<?) 



At the one-loop order, the trace of the vacuum polarization tensor is given 
by the fermion loop with effective propagators and vertices, as calculated in 
. The two-loop diagrams with their associated counterterms are shown in 
Fig. The justification of the use of counterterms is obvious. For example, 
when evaluating the first diagram in Fig. 1, one has to integrate over a region 
of phase-space where the gluon momentum L can be time-like (I? > 0) and 
hard (L^ ~ T). This contribution is already included when using effective 
fermion propagators in the one-loop diagram. The role of the self-energy 
counterterm is to remove from the two-loop contribution the part which is 
already included at one loop. A similar discussion can be made concerning 
the diagrams on the second line of Fig. |l[ which indicate that the vertex 
counterterm must be introduced. 



Figure 1: Two-loop contributions. A black dot denotes an effective propagator or vertex. 
Crosses are HTL counterterms. 




More formally, the counterterms automatically arise when we construct 
the perturbative expansion with effective propagators and vertices rather 
than with the original propagators and vertices. The effective perturbative 



2 



expansion is nothing but a reorganization of the usual perturbative expan- 
sion. The effective Lagrangian £ e fj, which leads to effective propagators and 
vertices, already includes some one-loop thermal contributions. In order to 
keep the same theory as that given by the original Lagrangian C one has 
to supplement the effective Lagrangian with counterterms to subtract order 
by order the loop-corrections already included in the effective Lagrangian so 
that 

C = C eS + C cX . (2) 

The counterterms are then treated perturbatively and thus appear in the 
expansion when looking at higher order topologies. 

The imaginary part of the retarded vacuum polarization in Eq. ([!]) can 
be expressed as a sum over different cuts through the diagrams of Fig. ||. 
Calculating the cut diagrams with the full complexity of effective propaga- 
tors and vertices would be a formidable task. To extract only the leading 
behavior of the two-loop contributions one can make some useful simplifica- 
tions. This leading contribution arises essentially when a hard momentum 
flows in the quark loop. This is a consequence of the phase space available, 
/ d^p ~ T 4 in the hard region, compared to g 4 T 4 for the soft region. Be- 
sides that, for soft fermion momentum P, the leading contribution should be 
contained entirely in the one-loop diagram since effective propagators and 
vertices give the complete behavior at scale gT. In the soft quark region, 
the two-loop diagrams (with associated counterterms) will only give sub- 
dominant terms while new dominant contributions can occur in the hard 
quark regime. Concentrating now on the region where the quarks are hard, 
it is sufficient to use bare vertices rather than effective ones and to keep 
only the time- like sector of cut fermion propagators, which largely domi- 
nates over the space-like sector. To see that, one can compare the behavior 
of the spectral density 5{P 2 )p = 0(l/p) in the time-like region with its be- 
havior g 2 T 2 J D ' jp in the space-like region when the momentum p is hard. We 
are therefore led to calculate the graphs of Fig. || with time-like cut fermion 
lines. 

With these simplifications, two types of cuts can be distinguished: those 
going through the gluon propagator and those going only through two fermion 
propagators. For the moment we ignore the latter because kinematical con- 
straints would require either P and R (see Fig. 2 for the notations) to be 
time-like and soft (which would lead to a non leading contribution because 
of the small size of the phase-space), or P to be space-like and hard so 
that the loop-corrections make these diagrams sub-leading compared to the 
one-loop result. Cutting through the gluon line, it is convenient to distin- 
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Figure 2: Two- loop contributions in the hard momenta limit. Lines without any dressing 
denote the hard limit of effective propagators, in which we keep a thermal mass. 

guish the space-like L 2 region from the time-like one. The case of L 2 < 
has been discussed in detail in || where it has been shown to contain a 
new physical process (bremsstrahlung) absent from the one- loop result. No 
counterterms were needed in this region. From now on, we concentrate on 
the case L 2 > 0. Physically, this region corresponds to Compton and anni- 
hilation processes, which are already included in an approximate way in the 
one-loop calculation of Braaten, Pisarski and Yuan (BPY) ||. We expect a 
better approximation of these processes to come out from the evaluation of 
the two-loop diagrams for the following reason: when evaluating these two- 
loop graphs, we do not neglect P and R compared to L in the self-energy 
and vertex corrections since P and R are hard and comparable to L. On 
the contrary, in the HTL approximation "external" momenta P and R are 
neglected compared to the loop-momentum L in the resummed self-energy 
corrections to the fermion propagators and in the estimate of the loop cor- 
rection leading to the effective vertex. We therefore anticipate that in our 
two- loop calculations, terms of leading order in L in the matrix element will 
be compensated by counterterms while terms of lower degree in L will lead 
to new contributions. 

The counterterm method described here is not the only way which can 
be used to calculate the two-loop diagrams while avoiding double counting 
of thermal corrections. Since, as mentioned above, effective propagators and 
vertices do not have the correct behavior in the hard limit (more precisely 
when some of their external momenta are hard and space-like) , an alternative 
way to proceed consists of introducing a cutoff scale intermediate between gT 
and T. In this method, one uses effective propagators and vertices in loops 
carrying a momentum below the cutoff. Above the cutoff, bare propagators 
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and vertices can be used, and a loop correction must be inserted on the 
hard propagator. When adding the soft and the hard contributions the 
cutoff dependence should drop out 0, This method has been pioneered 
by Braaten and Yuan in their calculation |l3|] of the rate of energy loss in a 
hot plasma, and used later in the calculation of the rate of hard real photon 
production by a quark- gluon plasma [14, 15]. However it was not followed 
in H where the effective propagators were used up to the hard momentum 
scale. 

A remark is worth making concerning the use of the cutoff method in 
[14, 15 1 . Indeed, in these two papers, a bare gluon propagator is used even if 
the cutoff does not constrain the gluon to be hard. As a consequence of the 
misuse of the cutoff in this way, [14, 15] missed the important contribution 
coming from bremsstrahlung J5J. 

In the following we will obtain analytically the leading \n(l/g) behavior 
of the vacuum polarization diagram. The argument of the logarithm is the 
ratio of a hard scale over a soft scale of order gT. Technically, this arises 
from integrals of type / dp/p or fdl/l. To extract this leading behavior 
one does not need the full details of the effective propagators in the soft 
region. It is sufficient to approximate the fermion and gluon dispersion 
relations by their asymptotic forms in the hard region: the error involves 
ratios of soft scales which are of order 1 compared to ln(l/<?). Furthermore, 
only the upper branches (quasi-particles) of the dispersion relations will 
contribute because the lower branches (collective modes) of both the fermion 
and the gluon decouple exponentially fast in the hard region and therefore 
cannot contribute a leading logarithm term. A last technical simplification 
consists in keeping a constant mass, independent of the momentum, along 
the dispersion curves: thus we have m? F ~ m 2 ~ g 2 T 2 , the exact expression 
of the thermal masses being irrelevant for the logarithmic behavior. 



3 One-loop result 

The virtual photon production rate has been calculated at one loop in the 
effective theory Q. In their calculation, BPY distinguished three types of 
terms depending on whether the pole (time-like) part or the cut (space-like) 
part of the fermionic spectral density is taken into account. We already 
extracted in |5| the leading logarithmic behavior of the cut-cut contribution 
and we showed that it should be compared to the L 2 < part of the two-loop 
diagrams. We give now the remaining leading logarithmic part associated 
to the pole-cut contribution (the pole-pole term does not lead to a loga- 
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rithm because hard time-like momenta are kinematically suppressed) which 
describes the Compton and the annihilation processes. From Eq. (11) of |J, 
taking the hard momentum limit in the integrand, one can easily extract 
the leading logarithm given by 



dN 



d 4 xdqod 3 q 



q2 (^)mHA) ■» 



Hoop f 



where N is the number of colors, C F = (N 2 - 1)/2N, m\ = g 2 C F T 2 /8 is 
the thermal mass of soft quarks and e/ is the electric charge of the quark of 
flavor /, in units of electron charge. 

For later use, it is convenient to translate this formula into an expression 
for the photon polarization tensor: 

Imn ^ (80 , o) ,_^^ (?e?)?0Tln (_^_). (4) 

4 The two-loop calculation 

4.1 General expression and the logarithmic behavior 

The two-loop expression (without counterterms) has been derived in ||] un- 
der the same simplifying assumptions as those used here (hard fermion mo- 
menta, no effective vertices). There is was found 

T tt ra u, \ NC f 2 2 M 4p f d ^ L r / \ / \ i 
Imn ^{q ,q) = — e g J J ^Kfro) -n F (p )\ 

x [n B (lo) + n F {r + l )]S(P 2 - m 2 F )5((R + L) 2 - ml)e( Po )e(r + l ) 
Traced i t Trace pff i 



x J2 Pa{L)P a pa {L) 

a=T,L 



self 



(5) 



R 2 (P + L) 2 R 2 R 2 
where e is the electric charge of the quark running in the loop and where 



the notation R = (ro, J r 2 + m 2 F r) includes the thermal mass shift on the 
fermion propagator: 



R 2 ± ir e R 2 - m 2 ± ir e ' 
After contracting over the transverse and longitudinal gluon projectors 



(6) 



13, 



p T m rrrr , (L p ~ l U p )(L a - l U a ) 

P pA L ) = 9pa ~ U p U a + j 2 (7) 
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(8) 



where U = (1, o) is the 4- velocity of the plasma in its rest frame we arrive 
at: 



£ Pa{L)P a pa {L) 
a=T,L 



Trace** 7 1 



Trace** 7 



+ 



(p T (L) - p L (L)) 



R 2 {P + L) 2 R?R 2 
4p 2 (cos 2 6>' - 1 



R 2 (P + L) 2 



L 2 — 2 



= Pa{L)\M a \ 
a=T,L 



{Q + Lf 



R 2 (P + Lf 
-2 P JL) (1-2 



(Q 2 p L (L) + L 2 p T (L) 



(Q • Lf 
r 2 (FTl) 2 



Q 2 L 2 



+ 



R 2 (P + L) 2 



+ 



(R 2 ) 2 {{P + L) 2 ) 2 



9) 



where 8' is the angle between r and I (or p and I since the virtual photon 
is static). Unlike in our previous work which dealt with the contribution of 
the space-like part of the gluon phase space, we consider here the time-like 
part and use: 

p a (L) = 2^e(/ )^ 2 -Ren a ). (10) 

Eq. (J5|) accounts for the cuts shown in Fig. 0. We still have to include 
the symmetric cut for the vertex diagram and also the diagram with the self- 
energy correction on the lower line. This will be done simply by multiplying 
the final result by an overall factor of 2. 

Since the external photon is massive, Q 2 ~ g 2 T 2 , there does not appear 
any collinear singularities of the type discussed in 0] and the logarithmic 
behavior arises only from terms like J dp/p or J dl/l. A simple power count- 
ing will then allow us to isolate the relevant terms in Eqs. (JsJ) and @j. We 
can use the following rules 

. Jd P0 5(P 2 -m 2 F )^-, 
• Jdl 5(L 2 -R e n a )~y, 



J d cos 9' 5((R + Lf 



m" ~ — 



1 

pi ' 



• n F (r ) - n F (p ) ~ — since q = r - p is soft, 
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n B (l Q ) + n F (r + Z ) ~ — ~ T , 



T T 
Zo ~ T 

to write, ignoring any irrelevant angular integrations, 

dl 



a=T,L 



Estimating i? 2 ~ 2pqo and likewise for (P + L)^, we find that only the terms 

, (Q ■ L f 

'R 2 (P + L) 2 
leading to 

Imn M / ~ e 2 g 2 q J dpdl (o(±) + 0(y)) (13) 

contribute to the leading logarithmic behavior. It has been checked by 
explicit analysis that all the other terms do not have the proper scaling 
behavior to give a logarithm. It can be noted here that the relevant terms 
for the present analysis are totally different from those in the bremsstrahlung 
case [1? < 0) . This is related to the different behavior of the gluon spectral 
density in the space-like region and the time-like region for hard momentum. 
More precisely, in our calculation of the trace using the Feynman gauge, the 
leading bremsstrahlung terms were all contained in the vertex corrections 
whereas for Compton and annihilation they appear in the diagrams with 
self-energy corrections. 

To summarize, in order to obtain the leading logarithmic terms for pos- 
itive L 2 it is sufficient to calculate 

d A P r d A L 
(2^7 (2^ 

x K(Z ) + n F (r + l )]5(P 2 - m 2 F )5((R + L) 2 - m 2 F )e( Po )e(r + l ) 

QL 



2p T (1 - 2 J 9 „ L) r 2 N , j -IV, \\.+Q.L\= I + 7==^ 1 j ■ (12) 




Imn W /(? ,o) ^4NC F e 2 g 2 J J [n F (r ] - n F (p )} 



x 2vre(/ )5(L -Ren T ) 




Elmn 1 ( go ,o)+Imn 2 (g ,o) (14) 

The above decomposition of Imll into two parts is natural: Imlli (the term 
1 in the square brackets) is dominated by hard p while the logarithm comes 
from the I integration; on the other hand, Im H2 (the terms proportional to 
Q ■ L) is dominated by hard I and it has a logarithm in the p integration, 
which indicates that we can make different approximations in each of these 
parts. From the discussion in section ^, we anticipate that Imlli is a new 
two-loop contribution while Imll2 should be compensated by counterterms. 
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4.2 Hard p region (Imlli) 

Being dominated by hard p, we can neglect m F in Imlli, and hence we can 
use the same kinematical approximations as in Q. Using the 6((R + L) 2 ) 
function, we get for the angle 9' between p and I: 

„, (r + l ) 2 -p 2 -l 2 
cos 9= . (15) 

Additionally, cos 9' must be kept within [— 1,+1], which places some con- 
straints on the phase space. Eq. (|l~5|), together with the condition on cos 9', 
gives the following inequalities: 

(l -l + Po + qo-p)(lo + l + Po + qo+p) <0 (16) 
(l -l+Po + qo+p)(lo + l + Po + qo-p) >0, (17) 

with \po | = p. We assume p fixed and hard and solve the inequalities for Iq 
and I which leads to the phase space reduction seen in Fig. [3|. The different 
allowed regions in Fig. ^ can be interpreted in terms of physical processes 
(Fig. H) determined by the relative signs of the energies po, Iq and ro + 
Iq. We have Compton processes in regions (I) and (III), quark-antiquark 
annihilation in region (II) , and plasmon decay in region (IV) . The symmetry 
of the integrand in Imlli with respect to the change of variables P «-> —R — 
L, indicates that regions (I) and (III) give equal contributions. Therefore, we 
will consider one region (I for instance) and multiply the end result by 2. The 
domain of integration of Eq. (|l4|) is constrained by the S function Eq. ( |l0|) , 
which further reduces the allowed phase space to the intersection of the 
dispersion curves with the regions I to IV. We also note that the longitudinal 
modes of the gluon cannot give a large logarithm in the momentum integral, 
since longitudinal modes are exponentially suppressed in the hard region. 
The transverse gluon polarization tensor introduced in Eq. ( |Io[) is given 

by: 

x + 1 



H t (Iq, I) = 3m 



x 2 x(l — x 2 ) , 
h — — : -In 



1 



(18) 



where x = Iq/1. For a time-like gluon, H T (lx,l) is a function of x varying 
slowly between m 2 /2 and 3m 2 /2. Therefore, at the logarithmic accuracy in 
the evaluation of Eq. (|14|), we can take II T ~ m 2 . Making this approxima- 
tion, we obtain the following expression for Imlli : 

T _ , , 8NC F e 2 g 2 f , , , , f , 
Imlli(go 3 o) w —-3 — q J dpn F {p) J dx 
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Figure 3: Allowed domains in the (lo,l) plane for po = ±p. The area shaded in dark 
gray is excluded by the delta functions. The region shaded in light gray is below the 
light-cone (dotted lines) . The thin solid curves are the transverse dispersion curves of the 
thermalized gluon. The vertical dotted line is the separation between e(po)e(ro + lo) = +1 
and e(po)e(ro + lo) = — 1- 



m g) \ n B V x ) + n F (~P + <?0 + lx)\ (19) 

where we have used n F (po + qo) — n F (po) « qon' F (p). In the normalization, 
we have taken into account the above mentioned factors of 2. The integrals 
can easily be done to logarithmic accuracy and we obtain: 

^n,(.,o)-^(E4),T t (-|^) > (20) 

where we have reintroduced the summation over the flavors running in the 
quark loop. The occurrence of m g indicates that this logarithm cannot be 
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Figure 4: Examples of physical processes included in the phase-space of of Fig. Region 
I: po < 0, ro + lo < 0: Compton scattering of an antiquark. Region II: po < 0, ro + lo > 0: 
qq annihilation. Region III: po > 0, ro + lo > 0: Compton scattering of a quark. Region 
IV: po > 0, ro + lo < 0: plasmon decay into qq and a photon. For the gluon, we use an 
approximate dispersion relation, with a constant thermal mass. The quarks being hard 
need not be effective. 



found at one-loop in the HTL scheme since only the thermal fermion mass 
appears at this level. As we will see later on, this is not compensated by the 
counterterms. We find that the annihilation process (region II) contributes 
to Im IIi at leading order e 2 g 2 qoT, but without a logarithm since I is con- 
strained to be hard in this region (see Fig. |3|). Similarly, the plasmon decay 
process (region IV) can be ignored at the logarithmic accuracy. 



4.3 Hard / region (Imn 2 ) 

In order to calculate Im II2 , we follow the same procedure as in the previous 
section but interchanging the roles of I and p. Since Im II2 is dominated by 
the hard /-region, we can neglect the gluon thermal mass. Again we start 



with the 5 function 5((R + L) 2 — m 2 F ), which provides us with the angle 9' 
via the relation: 



-m 2 F + 2r l + L 2 

cos#= £ . (21) 

2rl v ' 

As before, we must enforce — 1 < cos#' < 1, which implies the following 
set of inequalities: 



4 

(r - r + l - l)(r + r + l Q + /) < m 2 p . ( 23) 

We can further simplify these inequalities and write them as: 



{r -r + lo + l)(ro + r + l -l)>m z F (22) 

2 

F 



\f (P - I) 2 + m 2 F <lo + r or l + r < -^J{p-l) 2 + m 2 F (24) 
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Figure 5: Allowed domains in the (po,p) plane for lo = ±Z. The area shaded in dark gray 
is excluded by the delta function constraints. The areas shaded in light gray are below 
the light-cone (dotted lines) . The thin solid curves are the mass shells of the thermalized 
quarks. The vertical dotted line is the separation between e(lo)e(ro + lo) = +1 and 
e(l )e(r + lo) = -1. 

-yjip + 2 + K < l + r < ^{p + If + ml . (25) 

Keeping I hard, the above inequalities lead to a reduction of the allowed 
domain in the (po,p) plane (see Fig. ||) where the case 1$ = I and Iq = —I 
have to be distinguished. 

As before, the various regions admit a physical interpretation. Regions 
I and II represent Compton scattering on an antiquark and on a quark 
respectively, while region III contains quark-antiquark annihilation. Region 
IV describes quark decay into quark, gluon and photon, but it is not allowed 
if the initial and final quarks have the same masses. We note also the absence 
of plasmon decay, which is entirely due to the different approximations made 
in this section and in the previous one: here we assume L 2 = but keep the 
fermion mass and this forbids the decay of the gluon into a photon and a 
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massive quark-antiquark pair. 



III III 




Figure 6: Examples of physical processes included in the phase-space of Fig. ^ Region 
I: po < 0, ro + lo < 0,lo < 0: Compton scattering of an antiquark. Region II: po > 0, 
ro + lo > 0,lo < 0: Compton scattering of a quark. Region III: po < 0, ro + lo > 0,/o > 0: 
qq annihilation. Region IV: po < 0, ro + lo < 0,lo > 0: quark decay into qg and a photon. 
For the quarks, we use an approximate dispersion relation, with a constant thermal mass. 
The gluon being hard is treated as a bare gluon. 



The change of variables P «-> —R — L in the Q ■ L/(P + L) 2 term in 
Imll2 allows us to write: 

/cfiP f d^L 
J^Ts J 7^pK( r o)- n FW] 



x K (*o) + n F (r + l )]5(P 2 )5((R + L) 2 )e( Po )e(r + l Q ) 



x 2tt€(1 )5(L 2 ) 



2 .Q-L 



R 2 



(26) 



This integral can be evaluated relatively easily. Consider, for example, region 
III (Iq = I, po < 0, ro + Zo > 0, i.e. qq annihilation). For I fixed and hard, 
the fermion dispersion curve po 
of the region at 



p 2 + m? intersects the boundary 



Pn 



SI 



2q 



Pn 



I 



(27) 



so that, after doing all angular integrations and using the 5 functions, one 
finds 



ImIIn 



e 2 9 2 
(2vr) 3 



OO Pma::, 

Idl 



pdp 



o 



x[n B (l) + n F (l-uj p )] 



[n F (u) p ) - n F (u p - q )} 
2q uj p - ql 
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NC F eV 
2 (2vr) 3 



Idl I dtUp 



•> 2 

F 



290 

x[n.(0 + n,(l-^)l ^-'''^-»>' (28) 

We can neglect in front of / in the term n F (l — uj p ), and the two integrals 
decouple. The integral over I leads to the usual hard thermal factor 7r 2 T 2 /4 
while the integral over p yields a logarithmic factor ln(qoT / 'm F ) . Using the 
same method, one finds that the contribution of the sector Iq = —I in Fig. || 
is similar to Eq. (|28|) so that, to the required accuracy, Eq. (]26|) leads to 

i m n 2( ,„,o) « -^^(VJ^th, (29) 

where we have included the charge factor of the quarks. The normalization 
factor takes into account the contribution of all the ways of cutting through 
the photon polarization diagrams. 



5 Counterterms 



We have calculated the two- loop diagrams using various approximations. 
As a consequence, in order to be consistent, we have to use the same ap- 
proximation when evaluating the counterterms. In the appendix, we show 
that the complete computation of the counterterms gives the same logarith- 
mic contribution as the simplified version of the counterterm diagrams we 
discuss in this section. 

The simplified version of the vertex counterterm and the self-energy 
counterterm are depicted on the right of Fig. 0. 

Their general structure is as in Eq. (||) but the integrals have now to 
be evaluated using the HTL approximations, namely L hard with I? = 
and P and R neglected with respect to L with, as before, the notation 
P = (p , yjp 2 + m F p), R = (p + Qo, \Jp 2 + m F p). This implies (R + L) 2 
2R ■ L. One then writes 



Imi% ,o) 



NC 



ct 



F 2 2 

eg 



d A p 



d A L 



(2tt) 3 J (2tt) 



bvOo) - n F (p )] 



x [n B (l Q ) + n F (l M Po )8(Py(2R ■ L)5{L?) 
x 27rg pa 



Trace HTL 



Trace^Jseif 



(2P ■ L)R 2 



_|_ sn J. i 



R 2 R 2 



(30) 



14 



where the traces can be evaluated using the Feynman gauge for the gluon. 
In this equation we make use of ro + Iq ~ Iq so that some e-functions drop 
out. 

For the vertex we find simply in the HTL approximation 

5p(7 Traced | vcrtcx « 32 R.L P.L = 0. (31) 

Hence the vertex counterterm is vanishing within the approximation used. 
This is not a general result but a consequence of using the Feynman gauge. 
The self energy counterterm is non vanishing and we find 

5p(7 Traced | self = 16 (2(P -R)(R-L)- R 2 P ■ L) 

« 16 Tf Q-L . (32) 

In the last line, we have neglected terms which are suppressed by inverse 
powers of p at large p which cannot contribute to a logarithmic factor in 
the p integration. Plugging this expression in Eq. (pOj) and comparing to 
Eq. (p6|), we find that these two equations coincide at the logarithmic level, 
except for the sign. Indeed, their sum receives a contribution only when p 
is hard and therefore does not have a logarithmic behavior. One concludes 
that, to logarithmic accuracy, 

Im n 2 (g , o)+Im II( g o,o) U ^0 (33) 

so that, as expected, the counterterms cancel the two-loop contribution aris- 
ing from the hard L phase-space. 

6 The total two-loop contribution 

We summarize in this section the complete results of the calculation of the 
virtual photon rate up to two loops in the effective perturbative expansion. 
We have to collect the following pieces: the one-loop result of || given in 
Eq. ||, the two-loop bremsstrahlung contribution derived in j|] and finally the 
contribution of the Compton and annihilation processes we have calculated 
in the previous sections of the present paper. Adding all three contributions, 
we obtain to leading order in ln(l/g): 
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+? ln M^)+ 21n H^)i • (34) 

The term in rr? / q 2 . arises from the cut-cut term in the one-loop calculation 
while the term in m 2 /qQ describes the bremsstrahlung coming from the 
two-loop diagrams. The last line combines the pole-cut contribution at 
one-loop, Eq. as well as the L 2 > pieces at two-loop, Eqs. (|20|), 
(f29|), and the associated counterterms Eq. (|33|). In the above equation, the 
choice of hiding some powers of gT in thermal masses is somewhat arbitrary. 
Nevertheless, it should be noted that, even if the four terms have the same 
order of magnitude when qo is soft, two of them are proportional to e 2 g i and 
therefore would be obtained only at three loops in the naive perturbative 
expansion. 

If we consider the leading term in an expansion in powers of g at fixed 
qo (which is precisely what one is doing in the bare theory), we should be 
able to recover the results of two-loop bare calculations performed by ||. 
Taking this limit in Eq. (|34|), we find: 

dN a 2 2 x f NC F g 2 T 2 \ ,k 2 n T 

?^^»^(i:4)(^f-j( T +2)M-). (35) 

which is precisely the result found in 11. One may notice that, had we kept 
only the terms calculated by BPY || , the g — > expansion of Eq. (|34]) would 
not reproduce the complete bare result. 

One may also wonder what happens when considering the ultra-soft pho- 
ton limit, qo <C m F or go <C m g , for which the pole-cut contributions in 
Eq. ([M|) seems completely wrong. A glance at Fig. ^ and Fig. [| shows that, 
in that case, the intersection of the dispersion curves with the boundaries 
of the physical regions are pushed to hard / or p values so that no logarithm 
is generated. This is reflected in Eq. (|34[), where some logarithms become 
small and eventually negative when qo goes to zero, indicating that the terms 
we considered no longer exhibit a large logarithmic factor. More technically, 
in the case of Eq. (^) for example, the integration generates terms such as 

ln(1 - exp( -^^ )) 

which reduce to the usual logarithm for the generic case qo ~ m g ~ gT but 
which lead to an exponentially suppressed factor when q <C m 2 /T. 

For completeness, we briefly mention here the case of thermal production 



of hard or very hard real photon (qo > T) M, 15]. It was found in the 2- 



loop approximation M that, as in the case of soft virtual photon production, 
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leading contributions arose from the space-like part of the gluon spectral 
density [l? < 0) associated to the bremsstrahlung process as well as to 
q — q annihilation with scattering. To complete the calculation, the time- 
like {I? > 0) part should be considered. However it is not necessary to do 
the calculation using the counterterm method advocated above since in the 
published works [14, 15] the two-loop contribution has already been included 
using the cutoff method. The final result contains a factor ln(qoT/m 2 F ) 
similar to those of Eq. (F 



7 Conclusions 

In this work, we complete the study of soft static lepton pair production in 
a quark-gluon plasma, in the two-loop approximation of the effective per- 
turbative expansion. At one loop, the result shows a logarithmic sensitivity 
to hard momenta in the loop. Since the extrapolation of hard thermal loops 
is not accurate enough in the hard space-like region, we anticipate that the 
one-loop approximation may not yield the complete result at the leading 
order in the expansion in powers of the coupling constant. We should then 
consider two-loop diagrams, and an explicit calculation shows that, indeed, 
they give a contribution at leading order. Two types of terms can be distin- 
guished, which are associated with different physical processes: 

(i) two-loop leading corrections to processes already included at one-loop. 
They are the Compton and annihilation production of the virtual photon, 
which are calculated incompletely at the one-loop level. 

(ii) New processes, not contained in the one-loop approximation, namely the 
bremsstrahlung production of the photon. This was studied elsewhere [||. 

Dealing with two-loop diagrams in the HTL resummed theory requires 
some care. Two approaches are possible. One can use the cutoff method 
where the two-loop corrections are taken into account only above some cutoff 
value of the loop momentum. But this is not enough to evaluate the two- 
loop diagrams using bare propagators and vertices only. On the contrary, 
effective propagators and vertices are necessary when appropriate. Failing 
to do this, one could miss the "new" contributions referred to above. 

As an alternative to the cutoff method, one can construct the perturba- 
tive expansion from the effective Lagrangian, and counterterms must then 
be included to avoid double counting. These are crucial to correctly cal- 
culate the physical processes which appear at different orders of the loop 
expansion of the effective theory. This is the approach advocated in this 
work. 
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Since the new production mechanism of virtual photon at two-loop order 
shows a logarithmic sensitivity to hard space-like momenta in the loop, it 
cannot be claimed that our result is complete, for the same reason that the 
one-loop prediction could not be trusted. The contribution of three-loop 
diagrams should be considered. However, since in the present calculation 
there is no physical process that would arise at four loops or more in the 
bare theory, leading contributions coming from more than three loops are 
not expected. 

A Counterterms 
A.l Preliminaries 

In this appendix we give the exact evaluation of the counterterms. We show 
that the logarithmic behavior coincides with those found in the paper using 
the simplified version of the counterterm diagrams. The result according to 
which the HTL part of the effective vertices does not modify the calculation 
of the logarithmic part is to be expected on the basis of a quite general 
argument. Indeed, we know that the HTL correction to the qq'j vertex 
behaves like: 



where L = (1,1) and where P and R are the momenta of the quark and 
of the antiquark. This means that this HTL correction has a very different 
scaling behavior when P and R become hard, compared to the bare part of 
the effective vertex (which is independent of P and R). As a consequence, 
we do not expect a contribution to the logarithmic part from the HTL 
correction to the vertices. This is what we check explicitly in the following. 

To calculate exactly the counterterms of Fig. |l|, we follow H, where the 
authors relate the imaginary part of the effective vertex to that of the effec- 
tive propagator, in order to write the final result in terms of the propagator 
spectral density only. In this paragraph, we define some notations, and some 
useful formulae. 

The effective propagator can be written as: 



with D T (p,x) = ~[p 2 (x — t) — rr? tau + m 2 F (TX — l)Qo(x)] ,(37) 




(36) 




1 



P 
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where p T = 7 — T7 • p, and Qo(^) = \ hi((x + l)/(x — 1)) is the first Leg- 
endre function. The presence of Legendre functions in the vertex and the 
propagator allows us to relate the imaginary part of the vertex (or the imag- 
inary part of the product of a propagator with a vertex) to the imaginary 
part of the propagator. The spectral density of a quark and the imaginary 
part of the product of a vertex with a propagator are given by: 

I p 2 (x 2 — 1) 
Pt(p,x) = ~ 2Im j) fa x \ = — ^2 [S(px - uj±(p)) + 5{px + uj T (p))] 

+(3 ± (p,x)8(l-x 2 ) , (38) 



Qo(x) \ _ 1 (p\x-t) -m 2 F T 
D T (p,x) 



Im ) = 7- ^\ ' prM (39) 



with 

(3±(p,x 




x T 1 - (m F /p) 2 [Q (x) T Qi(a?)])^ + [(nm 2 F /2p 2 )(l T x)] 2 

(40) 

where Qi(x) = xQo( :c ) — 1 is the second Legendre function, and to±(p) are 
the solution of the dispersion equation D T (p, oj±{p)/p) = 0. 

We need also a formula giving the imaginary part of the product of two 
vertices and a propagator, given by: 

l m I Qo( x ) ) = P^Qo( x ) _l( p 2 {x-T)-m 2 F r \ 2 | . 

\D T (p,x) I ra 2 F irx — 1) ' .-■>,- 1 

where in the retarded prescription {x 

ImQo(x) = -'-0(1 -x 2 ) . (42) 

A. 2 Vertex counterterms 

There are two counterterms diagrams associated with the 2-loop vertex di- 
agram: one is shown in the lower right corner of Fig. [l] and the other one 
with the counterterm in place of the other loop. In what follows, we denote 
by —R the momenta circulating in the upper quark of the diagram of Fig. [l] 
and P to be the momenta of the lower quark. 

Applying the previous formulae, and following |Iq| , we get: 

o N +00 +00 +00 

dN ±N(j: f e 2 ) a 2 f r r 

' P dp dpo dr n F (p )n F {r ) 




d 4 xdqod 3 q 



3tt 4 



,2 
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x <%o - PO - r„) {2(1- P \^P) PAP)P-(R) 





2<?oP 






-2p 2 


- 2m 2 




2<?oP 




Po + r o 


-2p 2 


- 2m 2 




2<?oP 




P§ + »o 


-2p 2 


-2m 2 



2 



i 1 - - v j^v-^) 

2 

771 

+% 2 - pg)-^y(l " X 2 )[(l + + (1 - X)p_(i2)] 

" -o)^4(l - f?V)[(l + -)P + (P) + (1 - -)P -(HI } ( I-'!) 



8p<?o VP/ P P 

Since we want to look for possible double counting in Compton and 
annihilation processes, we must keep at least one time-like quark momentum 
(R 2 > or P 2 > 0). By symmetry, we can limit ourselves to one of these two 
regions, and multiply the end result by a factor of 2. Doing the logarithmic 
approximation (i.e taking (3±(p,x) m 2 /2p 3 (l =p x), ....), we deduce that 
the vertex counterterm vanishes at the level of approximation: 

dN 

~ , (44) 

vertex § 



d 4 xdqod 3 q 



which is equivalent to the result found using the simplified version of the 
counterterm diagram. 



A. 3 Propagator counterterms 

Making use of the same tools, we can calculate the counterterm diagram 
associated with the self-energy correction (diagram in the upper right corner 
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of Fig. |], with the other symmetric diagram where the counterterm insertion 
is on the lower quark propagator). 

Now, the relevant region of phase-space that we want to consider is the 
region where the line having a counterterm insertion is space-like while the 
other one is time-like. This is a consequence of limiting ourselves to the cut 
passing through the self energy insertion in the two loops diagram. 

After a lengthy but direct calculation we get (including the two possible 
diagrams of counterterms insertion): 



dN 



d A xdqod 3 q 



4iV(E /e 2) 



+00 +00 



of 



self 



3vr 4 



ql 



p 2 dp / dp / dr n F (p )n F (r ) 



x S(q -po- r )0(p 2 - p 2 ) 
W2II --'■■■^' r 



p 2 -rl 



+2 1 + 



Po~ r o 



p^d_ fp + (P) 
2 dp 



Zqop 



'Ydp 



MP) 



+ 1 + 



+ 1 



p 2 + r 2 - 2p 2 - 2m 2 F \ \ p 2 8 ( (3 + (P) 



iqop 

PO + r ~ 2 P 2 - 2m F 

Zqop 



2 dp \ p 

p^d_ fP-(P)\ 
'Ydp 



A_(P) 
A + (P) 

MP) 
J-MP) 



A + (R) 
A_(P) 



(1 - x 2 )[(l + x)A+(R) + (1 - x)A-(R)] 



P 



+2—x{x - 1) 1 
Qo \ 

-2—x{x + l) ( 1 + 



P 1 2 
— (x 

qo 



Zqop 

Po- r o 

Zqop 



MP)A_(P) 



P+ (P)A + (R) 



1)1 + 



Po + r o ~~ 2 P 2 ~~ 

Zpqo 



+ ^(x 2 -l) U-^^- 2p2 - 2m ' 



qo 



Zpqo 



MP)A+(P) 



/3_(P)A_(P) 



(45) 



where x = po/p and <9/<9p is the partial derivative with respect to p at 
constant x, and where A±(P) is the time-like part of the spectral density 
in Eq. (|39|), which can be approximated to logarithmic accuracy by 5(po =F 
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u + (p)) (i.e. the residue of the pole is approximated by 1, which is its value 
at large momentum). 

Finally, to logarithmic accuracy, we get: 



dN 



cfixdqocPq 



self 
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